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Spin-flip scattering of charge carriers in metals with magnetic defects leads to the low-temperature 
saturation of the decoherence time of electrons at a value comparable to their spin relaxation time 
r s . In two-dimensional (2D) conductors such a saturation can be lifted by an in-plane magnetic field 
By , which polarizes spins of scatterers without affecting the orbital motion of free carriers. Here, 
we show that in 2D conductors with substantially different values of the g factors of electrons ( g e ) 
and magnetic defects (gt), the decoherence time r v (B\\) (reflected by the curvature of magnetocon¬ 
ductance) displays an anomaly: It first gets shorter, decaying on the scale By ~ h/\gi — g e \fiBr s , 
before becoming longer at higher values of By. 

PACS numbers: 75.76.-t-j, 73.20.Fz, 73.43.Qt 


The electron interference results in a quantum correc¬ 
tion to the Drude conductivity and a positive magneto¬ 
conductivity (MC) of disordered metalsi In particular, 
the constructive interference of electron waves propagat¬ 
ing in time-reversed fashion along closed diffusive loops 
in two-dimensional (2D) conductors brings about a loga¬ 
rithmically divergent weak localization (WL) correction. 
In the absence of external magnetic flux piercing the elec¬ 
tron trajectories, this divergence of WL correction is cut 
off by the electron decoherence time t v . Application of 
the flux breaks the time-reversal symmetry, thus further 
diminishing the WL correction and leading to the low- 
temperature MC, cr(B±), where B± is the magnetic field 
component perpendicular to the plane of a 2D sample. 
The MC curvature, re = J^- 


oc rf, gives a mea- 


Bi=0 

sure for the electron coherence time in 2D conductors: 
doped semiconductor quantum wells, charge accumula¬ 
tion layers near semiconductor interfaces, thin metallic 
films, or atomically thin 2D crystals such as graphene 
and transition-metal dichalcogenides monolayers. 


The two leading decoherence processes at low temper¬ 
atures stem from the inelastic scattering of electrons off 
each other and off magnetic impurities. The electron- 
electron scattering results^ in the linear temperature de¬ 
pendence of the decoherence rate, 7y 1 = In 2 S/h ^ 

with —» 0. In contrast, the contact exchange inter¬ 

action with paramagnetic defects results in an apparent^ 
low-temperature saturation^— of the electron decoher¬ 
ence rate at t“ 1 (T —» 0) ~ r s _1 . The rate r” 1 charac¬ 
terizes the electron spin relaxation due to the spin flips 
in the course of electron scattering off randomly oriented 



FIG. 1. Dependence of the B± = 0 MC curvature, re oc on 
the in-plane magnetic field [see Eq. GDI. calculated for spin-i 
impurities with various g factors: red, yellow, and blue solid 
lines. Dotted and dashed blue lines demonstrate the magnetic 
field dependence of Cooperons without [Ci,i(By)-|-C'i,_i(By)] 
and with [— Co,o(By) + Ci,o(By) + Ci,i(0) + Ci,_i(0)[ spin- 
exchange processes for the case gi = 2. The inset shows the 
respective contributions to the full MC curve in the absence of 
an in-plane magnetic field [see Eq. m) and text for details. 
Here, T = 0.5 K, D = 100 cm 2 /s, and t s = 0.1ns. 


magnetic moments of impurities. 

It is common knowledge^— that electron spin relax¬ 
ation may be suppressed and extended by the polar¬ 
ization of magnetic impurities. In 2D conductors, this 
can be achieved by using an in-plane magnetic field, B y, 
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which polarizes the spins of the impurities but does not 
create any flux through the electron orbits. Measure¬ 
ments of k in various low-dimensional materials ^ 14 have 
shown a gradual increase of 77 (-By) associated with the 
spin polarization at gt^BBuS > T (here, gi and S are re¬ 
spectively the g factor and spin of a magnetic impurity; 
Hb = \e\h/2mc is the Bohr magneton). 

Here, we show that in some 2D materials the depen¬ 
dence of decoherence time on the in-plane magnetic field, 
t v (.B||), may be nonmonotonic: The naively expected 
polarization-induced increase of 77 , with B\\ is preceded 
by its decrease at weak fields (see Fig. |TJ) . This accelera¬ 
tion of decoherence comes from the precession dynamics 
of localized magnetic moments and requires the g factors 
of the impurities (gi) and electrons ( g e ) to differ from 
each other. To mention, if gi = g e , then the local mo¬ 
ments are static in the frame rotating together with the 
precessing electron spins, and in this case the 77 (By) de¬ 
pendence remains monotonic, being caused solely by the 
impurity spin polarization. For g t ^ g e , electrons witness 
the landscape of magnetic moments that varies in time 
with the frequency 


„ (ffe 9i)^BB || 

=- - -. ( 1 ) 

This temporal variation shortens 77 ,, if > 1. 

The latter condition is satisfied already at nonpolarizing 
fields, assuming that |1 — g e / gi\kTr s /hS 1. 

Polarization of the magnetic impurities renders them 
ineffective in the electron phase relaxation, thus lead¬ 
ing to a strong increase of the magnetoconductance sig- 
nal*i£~— Contrary to that, the effect of spin dynamics is 
quite subtle. We find the limitation on the magnitude 
of the corresponding change in the magnetoconductance, 
evaluate analytically the dependence of the 77 , on Bn , and 
relate it to the basic parameters of the itinerant electrons 
and magnetic impurities. 

To analyze the influence of spin-flip scattering in a 
conductor on the WL effect, we consider an electron 
wave propagating along a closed-loop trajectory, scatter¬ 
ing from disorder, V(r) = ]T r . U5(r - r*) + J2 r . J S, • 
cr<5(r — Tj). The Zeeman terms for electrons and impu¬ 
rities are —g e 9 BB\\a z and — gifisB^ (Sj) z , respectively 
(the z axis is chosen along the in-plane magnetic field), 
and a is the electron spin operator acting on the ±72 
spin states quantized along the z axis. 

To quantify the 77 (By) dependence, we express the 
WL correction 5a to conductivity^ in terms of two- 
electron propagators, “Cooperons” C aa i- tVV i(£, q;o/,w) 
[see Fig. Si)]: 


5a — —— [Co,o — Ci 0 — Ci 1 — Ci _i], 

Z7Trl 

Cs,M = Cs,M;CTCT , C'Cs,M;cro'') 


C = -"fVpT 3 / d£d 2 qn' F (E a ')C(£,q). (2) 



FIG. 2. Disorder perturbation theory diagrams for the WL 
correction to conductivity, (i) WL correction related to the 
Cooperon C arT i^i (e, q; ui',u>). Bold dots stand for the current 
operators and bold lines are the disorder-averaged Keldysh 
functions, (ii) Bcthe-Salpeter equation for the Cooperon. (iii) 
Combination of components of Keldysh functions involved in 
the kernel of the Cooperon equation. The first diagram in 
the right part of the equation contains all types of scattering 
described in Fig. Sa-d), while the second and third contain 
the spin-exchange processes (c,d) only. 


Here, Clebsch-Gordan coefficients = 

(5, M|Y 2 , <t; 1 / 2 , a') select from the Cooperon ma¬ 
trix C(e, q) = / duj C crcr '- 71J1 ’(£,q;u},uj) the singlet (5 = 0) 
and triplet (S = 1, M = —1,0,1) components defined 
in terms of the total spin carried by the two-electron 
propagator and its projection onto the external magnetic 
field B||. Also, 


£ a = £ - agifisB ||, (3) 

and n' F (e) = dnp(£)/d£ is a derivative of the Fermi dis¬ 
tribution function . 18 

The diagrammatic form of the Bethe-Salpeter equa¬ 
tion for the Cooperon matrix C is shown in Fig. m ii). 
Its important element is the disorder correlation func¬ 
tion represented by the dashed lines in Fig. [3] which is 
assumed to be short ranged and includes the following 
elements: 

(a) Correlator of spinless disorder, nuU 2 5(u})5(r — r'), 
where njj is the density of the pointlike potential scat- 
terers. 

(b) Correlator of the 0 -spin components of local magnetic 
moments that characterizes spin-dependent scattering of 
electrons without spin flip, nj J 2 (TKS z (t)S z (t'))5(r — r'). 
Here, nj is the density of the magnetic defects of spin S. 
The spin correlator (T KS z (t)S z (t')) is independent on 
the positions t and t' on the Keldysh contour, hence, it 
has only a Keldysh component with the Fourier trans¬ 
form njJ 2 (S 2 )5(uj)5(r — r'). 
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FIG. 3. Disorder correlation functions. First line: (a) 
Keldysh part of spin-independent scattering, (b) elastic spin- 
dependent scattering without spin-flip; (c,d) spin-flip scatter¬ 
ing, with the energy transfer in the spin-flip process shown 
along the dashed line. Second line: (c,d) Retarded/advanced 
correlator components. All (c,d) lines being implemented in 
Fig. [ 2 jiii) describe the spin exchange between two electrons 
and appear only in the Cooperons C 1 / 0 , 0 - 


(c,d) Spin correlators njJ 2 (T KS+(t)S-(t'))8{r — r'), 
where D(t,t') = (T K S + (t)S_(t r )) is mapped from the 
Keldysh time contour onto the matrix Keldysh space 
with components 

D r/a (uj ) = 2 i(S z )(uj - gni B B ± «0) _1 , 

D k (uj) = 4tt[S(S + 1) - (5 2 2 )]<5(u; - g^ B B). 

Here, (<!?") = [Z(a)]~ 1 d™Z(a), a = gijL B B\\/kT, and 

Z(a) = Xls e ° Sz the partition function for a para¬ 
magnetic scatterer. 

The thick solid lines in Fig. [2] stand for disorder- 
averaged electron Green’s functions G CT , obtained from 
the solution of the Dyson equation shown in Fig. [I] 

G R/A = (e a -v F ^ p +g e g B B\ l a±^(T- 1 + T~ 1 )j , 

Gf = [(1-2 n F (e a )\ ( G R - G A ), 

where r = lpTrjnuU 2 is the mean free time, t s = 
1/27T7 njJ 2 S(S + 1) is the spin relaxation time, 7 is the 
electron density of states, 

rp = [! - 2crMi (l - 2 n F {e- a ))] rp , 

M n = (Sp/S(S+ 1); f P «n F (|p| -p F ), 

and a = ± 1 /2 is the electron’s spin projection on the 
direction of the in-plane magnetic field. 

The spin structure of the Cooperons C\ ± 1 allows only 
for the (a) and (b) contributions to the dashed line in the 
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FIG. 4. First line: Diagrammatic representation of Dyson’s 
equation for the single electron Green’s functions calculated 
in the main order in U/vpft <C 1. Second line: Keldysh 
structure of the self-energy of the retarded/advanced Green’s 
functions. 


bottom row in Fig. [2jiii) , forbidding the spin exchange 
(c,d), and thus securing C aa '- aa ' (e, q; u/, w) cx <5(u/). 
The kernel of the Bethe-Salpeter equation for Cooperons 
Ci/ 0,0 includes spin-exchange contributions (c,d). Sum¬ 
ming up all three possible combinations of Keldysh func¬ 
tion components G A D K G R , G A D A G K 1 and G K D R G R 
in Fig. Oqii), where the frequency argument of the spin 
correlator is D R I K I A = D R l K l A [y> + gig B B), we get 

2 Y 1 ( D K + [1 - 2 n F (s- + u)]{D a - D r )) 

= 2wn J J 2 (s{S+ l)-<5 2 }-[l-2n F (e_)](5 z ))^(a;). 

The frequency dependence of this kernel enforces 
Coo'\-o,-a' (c q; a/, ui) oc 6(ui'), i.e. the energy trans¬ 
ferred through the impurity spin correlator can be only 
±giH B B ||, where the sign depends on whether the spin 
transferred to the defect is +1, or — 1 . After taking this 
into account, the equation for the Cooperon takes the 
form 


ihd T - r ) = 27r ^ r 2 ^( r ); 

Baa 1 — ^ari^a'rj 1 ~b Tr' )/2 T Tj* ^ T 
+ S a+(T I t r, +v i ^-4crcr , T s _1 M 2 + \<7 - T)\t~ 1 + 

+ i(cr - a')fl B ^j, (7) 

where D = v 2 f t/ 2 is the diffusion coefficient, C(e, r) is 
the Fourier transform of C(e, q) and Ax = (0,B±x) is 
the vector potential of the perpendicular magnetic field 
(note that B± -C B\\). 

Diagonalization of a matrix R produces Cooperons 
C'i.xi decoupled from all other Cooperon components 
and having decay rates 

ph = [l — M 2 =F Mi (l — 2 n F (e T ))] rp + Tp. (8) 

Here, e± [see Eq. ©] accounts for the energy transfer to 
an impurity in the process of spin-flip scattering. The 
Cooperon components Go,o and Gqo are coupled with 
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each other by spin-flip processes. The coupling generates 
combined modes decaying with the rates 


r o,± = [! + m 2 - M 1 (n F {£+) - n F (£-))] r. 


-l 


± 


i.+i 


U,-l - + Tj 


(9) 


Note that at By = 0 the average values (S z ) = 0 and 
(S*) = §5(5 + 1 ), so that t^Xi = Tq - = | T~ x + rf l 

and = 2t s -1 + rif , in agreement with earlier theo¬ 
ries 

Relaxation rates rrii and r 0 vary with By over 
two parametrically different field scales. For the 

scale is determined by the polarization of impurity spins. 
The polarization takes place at gifisB^S kT, and 
Eq. © then yields — rf 1 <C r s _1 . On the con¬ 

trary, the field dependence of is defined by the elec¬ 
tron and impurity spin dynamics. Under the condition 
|1 — g e /gi\kTT s /hS ^$> 1, the corresponding field scale is 
much smaller. Neglecting the spin polarization, we may 
simplify Eq. © to: 



As expected, the effect of the magnetic field depends on 
Us, the difference between the precession frequencies of 
the impurity and electron spins [see Eq. ©]■ The effect 
is absent if the corresponding g factors are identical. 

When substituted in Eqs. ©^ the four Cooperon 
modes obtained using Eq. © yield the WL correction 
to the conductivity at B± = 0 (the first term in square 
brackets comes from Ci,±i and the second from Ci/ 0 , 0 )) 


5a= 77 —p [ de Y, n' F (e a ) In + A a In , 
2nhJ a^± L r min T 0 ,+ _ 

A ± = ( T lX 1 - t t 1 )/( r o7+ - T o ,~) • 


( 11 ) 


Here, the ultraviolet cutoff under the logarithm is, typi¬ 
cally, T m i n ~ r, but for the description of the WL effect 
in graphene^ - — one should use for r min in Eq. (fill) the 
intervalley scattering time ri v , instead of the mean free 
path time r. The MC, studied as a function of B± for 
fixed B ||, takes the form 


cr(B ± ,B n ) -cr(0,H||) = -y-r f de n F (e a ) 

77 J a=± 


x { F 


B± 


B 


1,0;, 


F 




«<>,. 

1 1 


A Q }, ( 12 ) 


F(z) = \nz + if - + - , Bp >a = 


Ti/Ae 

DTp, a ' 


where if is the digamma function. 

One may see that the part stemming from Ci/ 0 0 
(square brackets) saturates at high B ±; it contributes 
(in units of e 2 /2 k h) at most (3/4) ln3 « 0.82 to the 


MC. This should be contrasted with the logarithmic 
growth with B± of the term coming from (the 

first term in braces). That peculiarity of the field depen¬ 
dence sets the dynamic range of MC useful for extract¬ 
ing the MC curvature, k(B ||) oc t 2 , using the expansion 
F(z) « z 2 /24 + 0[^ 3 ], as illustrated in the inset of Fig. |T| 
In the “high-temperature” limit, 

kT>hS/T s \l — g e /gi\, (13) 


the expression for n can be simplified further for all values 
of B ||, if inelastic e-e collisions are neglected (rf 1 —>• 0 ), 


(Dt „\ 2 ]2W b 4/3 I 

[ 9 M 2 (l + fn^) 2 /’ 

Wb = ch( 2 a) — ch(a), a = gi^sB\\/kT 1 


tt h \h/e 


Mi = 


(S + 5 ) coth[(5 + l)a] — § coth(la) 
5(5 + 1) : 


(14) 


where (S + 1)M\ is the Brillouin function [see Eq. ©]. 
The first term in braces in Eq. (11411 comes from Cooper- 
ons Ci,± 1 of Eq. © and has two asymptotes: 3 + 
29 +8 S(i+s) at B|| ^ kT/giUsS and 1(5 + 

l) 2 exp(2p,^sBy/AiT) at B\\ kT/g^s■ The latter ex¬ 
ponential asymptote is cut off by inelastic e-e scattering 
resulting in/c= if/rlli/e") ■ The second term in braces 
originates from C7i/o,o and its contribution to «;(By) de¬ 
cays with increasing By. Together, the two contributions 
provide the non-monotonic dependence of magnetocon¬ 
ductance curvature (which is conventionally considered 
as the measure of coherence time) over the in-plane field 
scale By ~ 7i /|gi — g e \l^BTs- This nonmonotonic depen¬ 
dence includes a local maximum at By =0 and a mini¬ 
mum at By ~ ^(|l-ff|^)” 2/3 , which is followed 
by the increase of k(B y) due to the polarization of defect 
spins [cf. Eqs. (fl3ll and (TT4l) ]. Eventually MC curvature 
saturates at the scale set by the inelastic electron-electron 
scattering decoherence time tt. 

The above-described anomalous behavior of the deco¬ 
herence rate occurs only when the electron or/and mag¬ 
netic defect have g factor values different from the free 
electron g = —2. The values of gi ^ —2 may be caused 
by the crystalline anisotropy effect on a heavy ion embed¬ 
ded in a 2D metal or semiconductor ( e.g ., graphene). For 
example, crystalline anisotropy splits states of a spin-| 
magnetic ion into two Kramers doublets with spin pro¬ 
jections ±2 and ±| onto the direction perpendicular to 
the plane of the 2D electron system. Then, the Zeeman 
splitting by the in-plane magnetic field realizes the cases 
of gi = —4 (yellow line in Fig. © and g z = 0 (orange line 
in Fig.© for the two doublets, respectively.— The case of 
g e = —2, < 7 i ^ —2 was realized in graphene exfoliated on a 
Si0 2 /Si substrate^ where the non-monotonic magneto¬ 
conductance is well-described^ by the theory presented 
here. Alternatively, the situation g e = 0 (green line in 
Fig.© can appear in p-doped transition metal dichalco- 
genides MoS 2 , MoSe 2 , WS 2 or WSe 2 , where, due to a 


















5 


large spin-orbit splitting, Kramers doublets of the hole 
states correspond to opposite spins in the opposite val¬ 
leys, and the external in-plane magnetic field does not 
lift the Kramers degeneracy^ 

In conclusion, the difference in the precession frequen¬ 
cies of the electron and impurity spins results in a non¬ 
monotonic dependence of the electron decoherence time 
on the magnetic field causing the precession. We find the 
magnitude and the functional form of that dependence, 
and relate it to the parameters of the itinerant electrons 
and magnetic impurities. Despite being small, the effect 
is important, as a manifestation of the very basic physics 


of magnetic moments in solids. 
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